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ABSTRACT 


This  paper  considers  a population  (composed,  for  example,  of  people, 
machines,  or  livestock)  subject  to  a randomly  occurring  defect  or  disease. 
If  there  exist  testing  procedures  capable  of  detecting  the  defect  before 
it  would  otherwise  become  known,  and  if  such  early  detection  provides  bene- 
fit, the  periodic  administration  of  such  a test  procedure  to  the  members 
of  the  population,  i.e.,  a mass  screening  program,  may  be  advisable.  An 
analytical  model  of  a mass  screening  program  is  developed  and  analyzed. 

In  the  development  of  the  model,  the  defect  arrival  process  is  Poisson 
and  the  ability  of  the  test  to  recognize  a defect  is  not,  in  general,  as- 
sumed perfect. 

The  model  focuses  on  the  time  since  incidence  of  a defect  as  the 
factor  which  determines  the  reliability  of  a test  administration  and  the 
value  of  a detection.  The  objective  is  to  minimize  an  arbitrary  increasing 
function  of  the  detection  delay  which  represents  the  disutility  from  the 
occurrence  of  the  defect  and  its  detection.  This  disutility  is  also  rep- 
resented as  a function  of  the  type  of  test  applied,  the  probability  of 
success  of  the  test,  the  testing  intervals  and  the  arrival  rates  of  the 
defects  over  Q subpopulations.  A relatively  simple  expression  is  derived 
for  the  objective  function  and  a comprehensive  mathematical  program  is 
presented.  The  case  of  perfect  test  reliability  is  then  considered  and, 
from  the  class  of  "cyclic"  schedules,  the  optimum  schedules  are  shown  to 
be  equally  spaced.  For  a polynomial  disutility  function,  properties  of 
the  optimal  schedules  are  presented.  Finally,  a method  for  determining 
the  disutility  function  from  experimental  data  is  suggested. 


A MODEL  FOR  OPTIMAL  MASS  SCREENING 
AND  THE  CASE  OF  PERFECT  TEST  RELIABILITY 

There  are  many  situations  where  a defect  can  occur  randomly  among 
the  members  of  a population--either  a population  of  human  beings,  in- 
animate objects,  or  perhaps  livestock--and  once  present,  exist  and  de- 
velop, at  least  for  a time,  without  any  manifest  symptoms.  If  the  early 
detection  of  such  a defect  provides  benefit,  it  may  be  worthwhile  to  em- 
ploy a test  capable  of  revealing  the  defect's  existence  in  its  earlier 
stages.  (A  defect,  disorder,  or  disease  will  generally  be  considered 
simply  as  a defect;  and  the  word  unit  or  individual  will  refer  to  a mem- 
ber of  the  population.) 

Of  course,  continuous  monitoring  would  provide  the  most  immediate 
such  revelation.  But  considerations  of  expense  and  practicality  will 
frequently  rule  out  continuous  monitoring  so  that  a schedule  of  periodic 
testing--a  screening  program--may  be  the  most  practical  means  of  achieving 
early  detection  of  the  defect.  In  general  terms,  the  question  then  be- 
comes one  of  how  best  to  trade  off  the  expense  of  testing  which  increases 
both  with  the  frequency  of  test  applications  and  with  the  cost  of  the 
type  of  test  used  against  the  benefits  to  be  achieved  from  detecting  the 
defect  in  an  earlier  stage  of  development. 

The  benefits  of  early  detection  often  depend  upon  the  application 
considered.  For  example,  in  a human  population  being  screened  for  some 
chronic  disease  (cancer,  diabetes,  glaucoma,  heart  disease,  etc.)  the 
benefits  of  early  detection  might  include  an  improved  probability  of  ul- 
timate cure,  diminished  time  period  of  disability,  discomfort,  and  loss 
of  earnings,  and  reduced  treatment  cost.  If  the  population  being  screened 
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consists  of  machines  engaged  in  some  kind  of  production,  the  benefits  of 
early  detection  might  include  a less  costly  ultimate  repair  and  a reduc- 
tion in  the  time  period  during  which  a faulty  product  is  being  unknowingly 
produced.  If  the  population  being  screened  consists  of  machines  held  in 
readiness  to  meet  some  emergency  situation,  an  early  detection  of  a defect 
would  reduce  the  time  the  machine  was  not  serving  its  protective  function. 
This  process  of  inspecting  a sizable  population  for  defects  is  called 
mass  screening. 

The  expense  of  testing  includes  easily  quantifiable  economic  costs 
such  as  those  of  the  labor  and  materials  needed  to  administer  the  testing. 
However,  there  can  also  be  other  important  cost  components  which  are  more 
difficult  to  quantify.  For  example,  in  the  case  of  a human  population 
subject  to  medical  screening,  the  cost  of  testing  includes  the  inconven- 
ience and  possible  discomfort  necessitated  by  the  test;  the  cost  of  false 
positives  which  entails  both  emotional  distress  and  the  need  to  do  unnec- 
essary follow-up  testing;  and  even  the  risk  of  physical  harm  to  the  testee, 
e.g.,  from  the  cumulative  effect  of  X-ray  exposure. 

The  rationale  for  constructing  a mathematical  model  of  mass  screening 
is  to  provide  a conceptual  framework  within  which  a mass  screening  program 
might  best  be  designed  and  its  worth  evaluated.  Such  a design  must  deter- 
mine which  kind  of  testing  technology  will  be  used.  Several  candidate 
technologies  may  be  available,  each  with  different  reliability  character- 
istics and  costs.  In  addition,  the  frequency  of  testing  must  be  decided. 

If  the  target  population  can  be  partitioned  into  subpopulations  according 
to  susceptibility  to  the  defect  (e.g.,  by  age,  family  background,  time 
since  last  overhaul  (for  a machine),  etc.),  then  the  best  allocation  of 
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the  testing  budget  among  the  subpopulations  must  be  determined.  Lastly, 
a decision  must  be  made,  for  the  population  and  type  of  defect  under  con- 
sideration, whether  a mass  screening  program  is  justified  at  all.  It  is 
felt  that  the  above  determinations  are  best  carried  out  within  the  concep- 
tual framework  of  a cogent  model  of  mass  screening. 

It  is  the  purpose  of  this  paper  to  present  a reasonably  comprehensive 
model  for  decision  making  with  regard  to  mass  screening.  The  model  util- 
izes a time-based  approach.  But  rather  than  seeking  to  analyze  or  minimize 
detection  delay,  as  in  some  of  the  literature  on  this  topic,  the  objective 
function  is  an  arbitrary  increasing  function  of  detection  delay.  (Detec- 
tion delay  is  the  time  between  the  incidence  of  the  defect  and  its  detec- 
tion, regardless  of  whether  that  detection  is  the  result  of  a screening 
test  or  of  the  defect  becoming  self-evident.)  The  reason  for  choosing  a 
general  function  is  that  the  disutility  experienced  upon  the  delayed  de- 
tection of  a defect  may  well  vary  in  a highly  nonlinear  way  with  the 
length  of  the  delay. 

Another  objective  sometimes  used  for  inspection  models  is  to  maximize 
the  lead  time  where  the  lead  time  is  the  time  difference  between  the  time 
of  detection  via  a screening  test  and  the  time  detection  would  otherwise 
have  occurred  had  not  a screening  program  been  in  existence.  However,  in 
the  case  in  which  the  test  being  used  is  perfectly  reliable  then  it  is 
possible  to  consider  the  lead-time  criterion  in  terms  of  the  following 
simple  function  of  detection  delay.  Suppose  that  at  the  (possibly  random) 
age  T (measured  from  the  time  of  incidence  of  the  defect)  the  defect,  in 
the  normal  course  of  its  development,  would  become  manifest  even  without 
a screening  test.  Then,  if  the  defect  is  detected  at  age  t,  the  lead  time 
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gained  is  (T  - t)+.  Therefore,  if  sup  T is  finite,  the  function  D(t)  * 
sup  T-  E(T  - t)+  may  be  used  as  a disutility  function  which,  when  minimized, 
will  maximize  the  expected  lead  time. 

Consequently,  the  results  presented  in  this  paper,  based  on  an  arbi- 
trary disutility  function,  generalize  as  well  as  extend  results  in  earlier 
work.  In  addition,  new  results  are  presented. 

In  the  next  section  a brief  review  of  the  literature  is  given.  Since 
the  work  on  inspection  models  is  very  large,  only  the  most  relevant  papers 
are  discussed.  The  interested  reader  may  refer  to  the  surveys  by  McCall 
[ 19651  and  Pierskalla  and  Voelker  [ 1976 ] for  a more  comprehensive  review. 

In  the  third  section  the  model  is  formally  stated  and  an  expression 
derived  for  the  expected  disutility  per  unit  time  incurred  under  a regime 
of  uniformly  spaced  test  applications.  In  section  4,  the  model  is  spe- 
cialized to  the  case  of  a perfect  test;  i.e.,  when  the  test  is  adminis- 
tered to  an  individual  with  the  defect,  the  defect  will  be  detected  with 
certainty.  Under  this  assumption,  a regime  of  uniform  test  intervals  is 
proved  optimal  within  a wider  class  of  "cyclic"  testing  schedules  for  a 
single  subpopulation.  The  problem  of  allocating  a screening  budget  among 
various  subpopulations  is  also  analyzed.  As  part  of  this  analysis,  the 
explicit  screening  schedule  (r^,  r^ r^),  where  r^  designates  the  op- 

timal testing  frequency  for  members  of  the  jth  subpopulation,  is  obtained 


in  terms  of  the  subpopulation  specific  incidence  rates  N^\j  and  the  budget 
constraint  for  the  disutility  function  D(t)  - atm.  For  D(* ) convex,  the 


above  solution  (when  m 


1)  is  shown  to  provide  a bound  on  the  ratios  r^/r 


r 


The  last  section  provides  a technique  to  estimate  the  shape  of  the 


disutility  function  from  empirical  data  in  the  case  of  a perfect  test. 
The  Appendix  contains  the  proofs  of  all  results. 
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Finally,  it  should  be  mentioned  that  the  population  (or  each  subpopu- 
lation, in  the  instances  where  a heterogeneous  population  is  considered) 
is  assumed  to  be  of  fixed  size,  N,  and  the  defect  to  arrive  according  to 
a stationary  Poisson  process  with  rate  NX  (N  could  be  a very  large  but 
finite  number).  It  may  be  somewhat  more  realistic  to  set  the  defect  ar- 
rival rate  proportional  to  the  number  of  defect-free  units,  rather  than 
to  the  total  number  of  units  in  the  population.  However,  it  is  also  as- 
sumed that  no  defect  can  remain  undetected  longer  than  T*,  even  without 
any  screening  tests  being  given.  Hence,  T*N X is  an  upper  bound  on  the 
expected  number  of  undetected  defects  in  the  population.  If  it  is  the 
case  that  once  a defect  is  detected,  the  afflicted  unit  is  replaced  in 
the  population  with  a healthy  unit,  then  T*NX  represents  a bound  on  the 
expected  difference  between  the  number  of  healthy  units  and  the  total 
number  of  units  in  the  population.  For  X small,  as  would  be  the  case 
for  a relatively  infrequently  occurring  disorder,  this  should  represent 
no  difficulty.  Consequently,  it  is  assumed  that  X is  small  relative  to 
N,  which  is  consistent  with  the  examples  of  potential  applications  which 
have  been  or  will  be  mentioned. 

2.  Literature  Review 

Some  of  the  early  papers  which  have  a bearing  on  time  dependent  models 
of  mass  screening  considered  a problem  formulated  as  follows:  Suppose  a 

single  unit  is  subject  to  failure  and  failure  may  be  detected  only  through 
an  inspection.  Each  inspection  has  an  associated  cost  and  there  is  another 
cost  proportional  to  the  duration  of  undetected  failure.  For  this  problem 
Derman  [ 1 96 1 ] found  a minimax  testing  schedule  when  the  failure  distribution 
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is  completely  unknown  and  there  Is  a finite  time  horizon.  The  reliability 
of  the  test  is  a constant,  p,  not  necessarily  equal  to  one. 

Roeloffs  [1963,  1967]  addresses  essentially  the  same  problem  with 
p * 1.  However,  the  decision  maker  knows  a single  percentile  for  the 
time-to- failure  distribution  of  the  unit. 

Barlow,  Hunter,  and  Proschan  [ 1963]  derive  explicit  rules  for  cal- 
culating the  optimal  inspection  schedule  for  a broad  class  of  failure 
distributions  (now  assumed  completely  known  to  the  decision  maker).  The 
test  is  assumed  perfectly  reliable,  i.e.,  p * 1. 

Keller  [1974],  again  assuming  p * 1,  developed  an  alternative  compu- 
tational technique  utilizing  the  calculus  of  variations  to  characterize 
the  optimal  inspection  schedule.  Again  there  is  a fixed  cost  for  each 
inspection.  And  there  is  a cost  which  is  a function,  not  necessarily 
linear,  of  the  detection  delay.  In  older  to  use  the  calculus  of  varia- 
tions, Keller  assumed  that  the  intervals  between  inspections  are  suffi- 
ciently small  (compared,  say,  to  the  mean  time- to- failure)  so  that  the 
inspection  schedule  may  be  expressed  as  a density  function  mapping  the 
unit's  age  into  an  associated  testing  frequency. 

Kirch  and  Klein  [1974] , whose  paper  is  addressed  explicitly  to  a 
mass  screening  application,  seek  an  inspection  schedule  which  will  mini- 
mize expected  detection  delay  subject  to  a constraint  on  the  expected 
number  of  examinations  an  individual  would  incur  over  a lifetime.  The 
test  is  assumed  perfectly  reliable.  Kuhn- Tucker  methods  are  utilized 
in  solving  for  the  optimal  screening  schedule.  The  point  of  view  adopted 
in  this  paper  is  similar  to  that  of  Kirch  and  Klein  and  Theorem  3 is  re- 
lated (although  it  is  somewhat  more  general)  to  their  main  result.  One 
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respect  tn  which  the  approach  here  differs  from  that  of  Kirch  and  Klein 
is  that  this  paper  considers  several  subpopulations  each  with  its  own 
characteristic  incidence  rate  for  the  disorder.  The  idea  is  then  to  al- 
locate optimally  a fixed  screening  budget  among  the  subpopulations.  Kirch 
and  Klein  instead  take  a longitudinal  view.  An  individual  through  his 
lifetime  is  subject  to  different  probabilities  of  incurring  the  defect 
and  a screening  schedule  is  optimized  subject  to  a constraint  on  the  ex- 
pected number  of  examinations  over  a lifetime. 

McCall  [1969]  considered  the  problem  of  scheduling  dental  examina- 
tions under  the  assumption  that  the  time  between  the  incidence  of  a cav- 
ity and  the  scheduled  dental  examination  controls  whether  the  cavity  re- 
sults in  a filling  or  an  extraction.  Cavities  are  assumed  to  occur  ac- 
cording to  a Poisson  process.  As  a generalization,  he  permits  the  time 
required  for  a cavity  to  become  beyond  repair  (by  a filling)  to  be  a 
random  variable. 

Lincoln  and  Weiss  [1964]  studied  the  statistical  characteristics  of 
detection  delay  under  the  assumption  that  the  times  of  examinations  form 
a renewal  process  and  that  the  probability  of  detecting  the  defect,  p(t), 
is  a function  of  the  defect's  age,  t.  They  derive  open  form  equations, 
similar  to  renewal  type  equations,  which  relate  the  density  functions  for 
the  following  entities:  the  probability  of  detection  at  a test  applica- 

tion (p C t ) ) , the  time  until  the  defect  becomes  potentially  detectable  and 
from  this  time  the  forward  recurrence  time  to  the  first  test,  the  proba- 
bility of  the  event  that  at  a particular  time  a test  occurs  and  all  prior 
tests  had  failed  to  detect  the  defect,  and  the  detection  delay.  For  the 
two  special  cases  where  p(t)  is  a constant  and  where  p(t)  is  exponential, 
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the  moments  for  the  detection  delay  are  derived  in  closed  form.  For  uni- 
form testing  intervals  (and  general  p(«)),  the  distribution  and  moments 
of  the  detection  delay  are  computed.  For  p(t)  = 1,  they  solve  for  that 
testing  schedule  which  maximizes  the  time  between  tests  subject  to  a con- 
straint on  the  performance  of  the  screening  program  relative  to  detection 
delay.  Two  such  constraints  are  considered.  The  first  bounds  the  proba- 
bility of  detection  delay  exceeding  some  threshold  T.  The  second  bounds 
the  mean  detection  delay. 

Theorem  1 in  the  next  section,  which  gives  the  expected  disutility 
in  terms  of  the  test  interval  and  test  type,  is  similar  to  the  objective 
function  given  by  Lincoln  and  Weiss  although  their  approach,  as  outlined 
above,  is  different. 

Several  recent  papers  have  studied  the  statistical  characteristics 
of  the  lead  time  provided  by  a screening  program- -either  a one-shot  screen 
of  the  population  or  a periodic  screening  program.  One  of  the  earliest 
such  papers  was  that  of  Hutchison  and  Shapiro  [1968] . 

Zelen  and  Feinleib  [1969]  point  out  an  important  use  for  the  expected 
lead  time  in  comparing  survival  time  data  (measured  from  the  time  of  detec- 
tion) gathered  from  a screened  population  versus  that  gathered  from  a con- 
trol population.  The  expected  lead  time  must  be  subtracted  from  the  ex- 
pected survival  time  of  the  screened  population  before  a valid  comparison 
can  be  made  against  the  control  group's  expected  survival  time.  In  a later 
paper,  Zelen  [ 19 71]  derived  expressions  for  mean  lead  time  in  the  case  of 
two  successive  screens . 

Prorok  [1973]  has  established  properties  for  the  lead  time  in  the 
case  of  n successive  screening  tests  conducted  at  uniform  intervals.  He 
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also  derives  the  proportion  of  preclinical  cases  detected,  both  at  a given 
test  and  over  all  n tests.  The  test  is  assumed  to  have  a constant  relia- 
bility not  necessarily  equal  to  one.  In  addition,  he  derives  an  estimator 
for  the  mean  lead  time. 


3 . A Model  for  Mass  Screening 

Although  most  of  the  definitions  are  stated  prior  to  their  use  in 
each  section,  listed  below  are  some  notation  and  conventions  used  at 
various  times  throughout  the  paper. 


(i)  1 (t) 

m 

r 


1,m  m£l, (t) 

l7’  r J 


if  - * t s 
r 

otherwise 


m+1 

r 


n 

. (ii)  II  x » 1 for  n < 1 

i=l 

n 

(iii)  £ x ■ 0 for  n < 1 

i=l 


(iv)  [x]  is  the  largest  integer  not  exceeding  x. 

(v)  The  phrase  "increasing  function"  shall  mean  a strictly  increasing 
function. 


As  mentioned  previously,  the  objective  function  will  not  be  to  maxi- 
mize expected  utility,  but  rather  to  minimize  expected  disutility.  The 
disutility  associated  with  a particular  detection  will  depend  only  on  the 
elapsed  time  since  the  defect's  incidence.  The  notation  D(t)  will  express 
the  disutility  incurred  if  detection  occurs  t units  of  time  after  incidence 
D(‘ ) is  assumed  throughout  this  paper  to  be  a nonnegative  increasing  func- 


tion. 


In  those  cases  where  sup  D(s)  < «°,  there  is  a natural  relation  be- 

s 2 0 

tween  the  notion  of  a utility  function  and  D(* ) ; namely, 


U(t)  = sup  D(s)  - D(t)  . 
si  0 


U(t)  is  a decreasing  function  expressing  the  disutility  avoided  by  de- 
tecting a defect  t units  of  time  after  its  incidence. 

Initially,  only  a single  susceptibility  class  of  size  N will  he  con- 
sidered. The  results  will  then  be  generalized  to  an  arbitrary  number  of 
subclasses.  The  times  of  incidence  for  the  defect  in  the  population  are 
assumed  to  form  a Poisson  process  with  parameter  N\  and  is  designated  by 
the  sequence  {S^} , k = 1,  2,...  This  assumption  is  motivated  by  the  fact 
that  any  arrival  process  with  the  following  characteristics  is  a Poisson 
process:  at  the  time  of  an  arrival  there  is  almost  surely  (i.e.,  with 

probability  one)  only  one  arrival,  that  the  number  of  arrivals  in  a time 
interval  does  not  depend  on  past  arrivals,  and  that  the  number  of  arrivals 
in  intervals  of  equal  length  are  identically  distributed  (Cinlar  [ 1975] ) . 
For  many  applications  such  as  the  arrival  of  diseases  like  cancer,  heart 
trouble,  etc.,  this  assumption  is  quite  reasonable. 

Were  a screening  test  of  type  i to  be  administered  to  the  individual 
with  the  kth  defect  (i.e.,  kth  in  the  order  of  defect  incidence)  at  time 
S + t,  then  the  test  outcome  random  variable  is: 


Definition: 

YJ(t) 


if  kth  defect  is  not  detected  at  time  t 
after  incidence 

if  kth  defect  is  detected  at  time  t after 
incidence . 


(3.1) 
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k 

Since  a defect  cannot  be  detected  before  its  incidence,  Y^(t)  = 0 for 

k 

t < 0.  Notice  that  the  argument  of  Y^(* ) refers  to  time  relative  to  the 
incidence  of  the  defect.  If  the  population  is  screened  at  time  s,  where 
s is  absolute  time,  then  the  kth  defect  will  be  detected  if  and  only  if 
Sk  * s and  Yk(s  - Sk)  = 1. 

It  is  assumed  Y^(t),  y|(s)  are  independent  except  when  k = j and 
t = s,  and  that  P^Yk(t)  = 1^  depends  only  on  l and  t.  This  latter  assump- 
tion makes  possible  the 

Definition:  p£(t)  = p(yJ(0  = l)  . (3.2) 

The  function  p^(* ) describes  the  reliability  characteristics  of  the 
type  JL  screening  test,  i.e.,  how  likely  such  a test  is  to  detect  a defect 
as  a function  of  the  defect's  age.  Note  that  p^(t)  = 0 for  t < 0. 

In  this  section  the  screening  test  is  assumed  to  be  administered  to 
the  entire  population  at  the  times  1/r,  2/r,  3/r,...  The  testing  fre- 
quency r is  a control  variable.  In  the  next  section,  upon  assuming  per- 
fect test  reliability,  it  is  shown  that  the  above  schedule  of  uniform 

testing  intervals  is  optimal  within  a wider  class  of  "cyclic"  schedules. 

“lc 

The  random  variable  S . denotes  the  time  at  which  the  kth  defect  is 

r,Z 

k 

detected.  S . depends  on  the  arrival  time  of  the  defect  (S  ),  the  type 
r , L 

of  test  used  (j i) , and  the  testing  frequency  (r) . 

Definition:  Sk  ^ = min  |n/r )Yk^  - Sk  ^ = Ij-  (3.3) 

n=  1,2,... 

Given  the  application  of  test  type  l at  the  times  {l/r,  2/r,...}, 
the  disutility  incurred  by  the  kth  defect  is  D^Sk  £-Sk^.  To  represent 
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the  total  disutility  incurred  due  to  those  defects  which  occurred  in  the 
interval  [j/r,  (j  + l)/r),  it  is  useful  to  define: 


Definition: 


r.ji.J 


(Sk)  . 


The  following  theorem  provides  an  expression  for  E[B  . .]  in  terms 

r f*1  > J 

only  of  the  disutility  due  to  detection  delay  (as  expressed  by  D(»))  and 

of  the  reliability  of  test  type  l (as  expressed  by  p (•)).  In  addition, 

the  theorem  shows  that  E[B  . .]  = E[B  ] for  j = 0,  1,  2,... 

r , j r , * , u 


Theorem  1 


00  . r n-1 

£[Br  | J ’ J D(u)p^(u)  IT  [l  - p^(u  )]du 

* n-1  , m-1 


n- 1 
r 


for  j - 0,  1,  2, . . . 


This  theorem  yields  a relatively  simple  expression  for  the  expected 
total  disutility  in  any  interval  of  length  1/r  and  using  test  type  l with 
probability  of  detection  p^(*)*  This  expectation  is  used  in  the  objective 
function  of  a mathematical  program  to  determine  the  optimal  testing  fre- 
quency for  a mass  screening  program  for  a heterogeneous  population.  In 
order  to  develop  this  mathematical  program,  the  following  definitions  are 
useful . 

_ n-1 

Definition:  B . - lim  — 72  eTB  . .1  . 

r>*  n j-o  T’Z>i 

B is  the  long-run  expected  disutility  per  unit  time.  (The  factor  r 

T,l 

enters  the  definition  to  convert  disutility  per  unit  testing-interval  into 
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4 


\ 


per  unit  time.)  By  Theorem  1,  EfB^  ^ ■ EfB^  ^ Q]  for  j * 0,  1,  2,... 

Therefore, 


r,i 


rE(Br,».o) 


n«l 


n-1 


® . — n-1 

rN\  S r D(u)p . (u)  II  [ 1 - p . (u  - - )]du  . (3.4) 

J 1 m-1  1 r 


Notice  that  if  test  type  i provides  perfect  reliability,  i.e., 
p^(t)  = 1 for  t s 0,  then  (3.4)  gives 

1 

®r  l * rNX  j r D(u)du  • (3.5) 

0 


Now  consider  the  problem  of  selecting  testing  frequencies  and  test- 
types  for  each  of  Q different  susceptibility  classes  which  together  com- 
prise the  whole  population.  These  classes  may  differ  from  one  another  in 
the  number  of  units  they  contain,  in  their  defect  incidence  intensity,  and 
in  the  cost  per  test  application  to  an  individual  for  a particular  type  of 
test.  The  subpopulations,  however,  are  assumed  to  share  a common  D(* ) 
function. 

The  solution  to  (3.4)  for  the  expected  long-run  disutility  per  unit 
time  in  the  context  of  the  prior  portion  of  this  section  may  be  regarded 
as  applying  to  subpopulation  j where  j - 1,  2,...,  Q by  writing 


Bj,r(jM(j) 


r(J)N  X(j) 


Z 


i-1 


Ir<J>  ^(u)P£(j) <u)  n [1  - Pfcj)1" * 7nr)]da  - 

i-l 

r(j) 


This  permits  a multi-subpopulation  screening  problem  subject  to  a 


budget  constraint  to  be  formulated 

as  follows: 

Minimize 

(r(l) r(QM(l),...,jl(Q)) 

$ g 

^ J.r(JM(J) 

(3.6) 

such  that 

^ N c ,r  £ b 

jml  j .M(j)  J 

(3.7) 

r . > 0 j - 1 , . . . , Q 

(3.8) 

K j)  j = 1,...,  Q 

(3.9) 

where 


b 


cost  per  application  of  a test  of  type  £(j)  to  an 
individual  of  subpopulation  j , 

number  of  units  (or  individuals)  in  subpopulation  j, 
budget  per  unit  time. 


In  order  to  make  this  mathematical  program  even  more  comprehensive, 
it  is  possible  to  add  constraints  on  the  amount  of  testing  labor  available 
and  on  the  capacity  of  the  testing  facilities  in  terms  of  the  number  of 
arrivals,  the  frequency  of  testing  and  the  type  of  tests  used.  In  addi- 
tion, the  cost  of  false  positives  can  be  included  as  a part  of  the  test 
costs  in  inequality  (3.7). 

For  example,  a constraint  on  the  total  labor  available  is: 


£ 


V:M(J)rj 


s l. 


(3.10) 


and  on  the  total  testing  facilities  available  is: 


-15- 


for  each  type  of  test  A(» ) used  over  the  subpopulations  being  tested, 
where 

■■  amount  of  labor  needed  to  administer  test  type  4(j) 
to  an  individual  of  subpopulation  j , 

= amount  of  testing  facility  time  needed  to  administer 
test  type  4(j)  to  an  individual  of  subpopulation  j, 

■ total  amount  of  labor  available  to  administer  test 
type  j i(«)  per  unit  time, 

* total  amount  of  facility  time  available  to  administer 
test  type  £(•)  per  unit  time. 

4.  The  Case  of  Perfect  Test  Reliability 

In  this  section,  it  is  assumed  that  p^(t)  = 1 for  all  t a 0;  that  is, 

^ there  is  perfect  test  reliability.  If  a defect  is  present  at  the  time  the 

test  is  administered,  that  defect  will  be  detected  with  certainty. 

Unlike  the  previous  section,  uniform  testing  intervals  are  not  assumed. 
Rather,  such  uniform  intervals  (within  a subpopulation)  will  be  proven  op- 
timal within  a larger  class  of  "cyclic"  testing  schedules.  Once  this  is 
done,  the  problem  of  allocating  a screening  budget  for  a target  population 
among  Q subpopulations  which  comprise  the  target  population  is  analyzed. 

That  part  of  the  budget  allocated  to  a particular  subpopulation  controls 
how  frequently  the  screening  test  can  be  administered  to  that  subpopula- 
tion. This  budget  allocation,  in  turn,  is  determined  by  the  relative  in- 
cidence rates  of  the  disorder  and  the  relative  costs  of  testing  which 
r 

characterize  the  different  subpopulations. 

Until  now,  r was  used  to  designate  the  testing  frequency  and  1/r  the 
interval  between  tests.  In  order  to  consider  more  general  testing  schedules, 
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a different  notation  Is  needed.  Suppose  that  the  jth  subpopulatlon  is 
tested  according  to  an  m(j)-cycle.  To  be  specific,  the  intervals  between 
tests  are  successively  x^,  x^,...,  Letting 

j V j 

y^  - Lt  xJ  , the  successive  times  of  test  administration  for  subpopula- 
v-1 

tlon  j are  yj,  y{ ^(J),  jJ(j)  + y{.... 

Let  D^x|,...,  x^j be  the  long-run  expected  disutility  per  unit 
time  for  the  jth  subpopulation  under  the  above  testing  schedule.  The 
problem  then  becomes  the  following: 


Q -/  1 

Minimize  12  D(x", 

(m(l) .... ,m(Q) ) j-1  V L 

(x1  x*  ) 

^ 1"  ' ' ’ m(l); 


•’  Xi(j)) 


(xQ  xQ  ) 

'‘V*  * m(Q)' 


(4.1) 


such  that  12  N c m(j)/y  , . £ b 
j-1  J J ni(j) 


x|  > 0 ; t ■ 1 m(J);  j * 1 Q 


where 


Nj  - the  number  of  individuals  or  units  in  the  jth  subpopulation, 

Cj  » cost  per  test  application  per  individual  in  the  jth  subpopu- 
lation, 

b ■ bound  on  the  average  cost  per  unit  time  for  all  tests  ad- 
ministered under  the  screening  program. 

For  the  moment  consider  a homogeneous  population,  i.e.,  let  Q - 1, 
and  suppress  the  index  j in  the  above  cyclic  testing  schedule.  Under  the 
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r 


assumption  of  a perfectly  reliable  test,  if  the  kth  defect  occurs  before 

time  y , it  will  be  detected  at  test  time  y,  if  and  only  if  y , < Sk  s y . 
m i J i-1 

Therefore,  the  total  disutility  incurred  due  to  all  defects  arising 
in  the  interval  (0,  y ] is 


m 

s 

i-1  k-1 


2 D(y,  -sk)i 


(yi-i’yi] 


(sk) 


where  yQ  - 0,  although  a test  is  not  made  at  time  0. 

Similarly,  the  total  disutility  due  to  defects  arising  in  the  inter- 
val (vy , vy+y)  for  v - 0,  1,  2,...  is 
m mm 


m °o 

2 E D(vy  + y,  - Sk)l . , (Sk)  . 

i-1  k-1  m 1 ^m^i-l^V7!1 


(4.2) 


Taking  expectations  and  using  Lemma  A given  in  the  Appendix,  (4.1)  becomes 

vy  +y 

m » ■'m  Jt 

E XN  D(vy  + y - s)ds 

i-i  i m 1 

VVyi-l 

<n  pyi  m -Xi 

- E XN  D(y  - s )ds  - E XN  | D(t)dt  . 

‘-1  y 1 1-1  Jo 

7i-l  U 


Hence,  the  expected  disutility  due  to  defects  arising  in  the  interval 

[vy_»  vy_  + y)  ^or  V - 0,  l,  2,...  does  not  depend  on  v and 
m m m r > 


m Xi 

D(x  , ....  x ) - ^ E [ D(t)dt  . 
1 m ym  i-1  J 


(4.3) 


some 


It 


i, 


will  now  be  shown  that  for  any  schedule  such  that  x^  + x.  for 

j€{l,  2,...,  m},  the  alternate  schedule  [x  , x. x } where 

12  m 
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U1 

x,  - x.  - . . . - x - £ x./m  is  feasible  with  respect  to  the  program  and 
12  “ j-1  J 

provides  at  least  as  small  a value  for  the  objective  function.  Once  this 
is  accomplished,  among  all  "cyclic"  schedules  all  future  consideration  may 
be  restricted  to  testing  schedules  characterized  by  equal  intervals  between 
all  test  administrations  (within  a single  subpopulation).  The  following 
lemma  will  be  required: 

Lemma  1:  For  D(* ) an  increasing  function,  h(x)  * J D(t)dt  is  a (strictly) 

0 

convex  function  in  x. 

/ m v m 

Since  h(*)  is  convex,  hCx^)  - h^  £ x /mj  < £ h(x  )/m  for 

' v-1  v-1  v 

i - 1,...,  m.  Adding  these  m equations  together, 


m m 

£ h(x  ) < £ h(x  ) 
i-1  1 v-1  V 


r E < r s h(xt) 

^m  i-1  1 ym  i-1 


m m 

where  y - £ x.  - £ x,  . Using  (4.3) 
" i-1  1 i-l  1 


D(x^  , X2 , • • • » ^ D(x^  , X2  9 • • • t 


Hence,  the  schedule  (x, , . . . , x ) yields  lower  long-run  expected  disutility 

1 tn 

per  unit  time  than  does  (x, , . . . , x ).  Since  both  schedules  require  the 

l m 

same  number  of  tests  per  unit  time  (m/y  ) and  since  x.  - x~  - ...  - x , 

m 1 z m 

then  any  m-cyclic  testing  schedule  may  be  replaced  with  a uniform  testing 
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schedule  without  increasing  average  testing  frequency  or  expected  dis- 
utility per  unit  time. 

To  return  to  the  general  case  of  n subpopulations,  replace  the  sched- 
ule ^x| , . . . , with  a uniform  schedule  utilizing  test  intervals  of 

y^,../m(j)  for  j - 1,...,  n.  Since  none  of  the  subpopulations  use  more 

ml  J ; 

of  the  testing  budget  than  before  (the  average  testing  frequency  remaining 
the  same  for  each  subpopulation),  the  new  testing  schedule  is  feasible  if 
the  old  one  is. 

The  above  result  is  given  formally  by 

Theorem  2:  If  p(t)  - 1 for  t 2:  0,  schedules  with  equally  spaced  test 

intervals  are  optimal  within  the  class  of  m-periodic  schedules  for  m * 

1,  2,... 


Now  let  x^  - x|  - x^ 


. * x^ , then  by  (4.3) 
m 


5 ofx- 


J-l 


xj 

Q X N mO)  Xi 

t -J-J-  S D(t)dt 

j-l  ym(j)  i-1  n 


. £ f D(tM, 

j-l  m(j)xJ  i-1 


, Z f 
J-l  *j  i 


D(t)dt  . 


(4.4) 


With  uniform  testing  intervals,  it  becomes  possible  to  define  the 
testing  frequency  for  subpopulation  j as  r . * 1/x  . Note  that  (4.4)  is 
consistent  with  (3.5)  of  section  3 where  the  assumption  of  uniform 
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testing  intervals  was  made  and  the  problem  of  imperfect  test  reliability 
analyzed. 

By  (4.4),  the  mathematical  program  (4.1)  can  be  rewritten 


Minimize 
(^1 » r2 » • • • » ) 


Q r r 

S N \ r J D(t)dt 

3-i  J j J ; 


such  that 


Ti  N c r s b 

^ J J J 


and  r j > 0 for  j - 1,  2 Q J 

where  the  constraint  equation  was  transformed  using 

m(J)/ym(j)  ■ m(j)/^x|m(j)^  = 1/xJ  » l/xj  = r . 

The  solution  of  (4.5)  will  optimally  allocate  the  screening  budget 
among  the  Q subpopulations  so  as  to  minimize  the  long-run  expected  dis- 
utility per  unit  tine  over  the  entire  population.  (If  r must  be  integer 
valued,  then  just  add  this  restriction  and  solve  (4.5)  as  a constrained 
knapsack  problem.) 

For  the  continuous  case,  let 


r J'dO 


for  r > 0 


r J D(t)dt  for  r * 0 . 


At  this  point,  assume  that  D(* ) is  a differentiable  function  on 
[0,  «).  Therefore,  f (• ) is  continuously  differentiable  on  (0,  ®)  and 
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f,(r)  " df  " d7[r  JrD<s)ds]  “ JrD(s)ds  -^D(±) 


, iD(Y)  .ID(I)  < o , y€(0,  i),  r > 0 


(4.6) 


where  the  last  equality  follows  from  the  mean  value  theorem  for  integrals. 
Also,  since  D' (•  ) > 0 


-£UtD<^-TD<T>] 

0 


“3  D'  (^)  > 0 


for  r > 0 . 


(4.7) 


So,  on  the  interval  (0,  «),  f(* ) is  a convex  decreasing  function  in  r. 

The  same  may  therefore  be  said  of  the  objective  function  in  (4.5).  The 
fact  that  the  objective  function  is  decreasing  in  each  of  its  components 

implies  that  the  constraint  £ N c r £ b is  binding.  Consequently,  this 

j-1  J J J 

inequality  may  be  replaced  with  equality  in  (4.5). 


Lemma  2:  If  lim  D(s)  - »,  then  lim  f(r)  =■  ® . 

s - ® r~  0 


Theorem  3:  If  lim  D(t)  * <*>,  then  there  is  a unique  solution 

t CO 

(r*,  r* r*)  to  (4.5)  which  satisfies  Xif'(r*)/ci  » Xjf'(r*)/Cj  for 

i,  j - 1,  2,...,  Q. 

In  the  proof  of  Theorem  3 in  the  Appendix,  it  can  be  seen  that  the 

hypothesis  that  lim  D(t)  ■ <*>  was  needed  only  to  establish  (A. 3).  There- 

t-« 

fore,  an  alternate  theorem  can  be  given. 

Theorem  31  : If  (r*,  r*,...,  r*)  is  a unique  solution  to  (4.5)  and  if  r*>  0 

for  J - 1,...,  n;  then  X1f'(r*)/ci  - Xjf'(r*)/Cj  for  i,  j - 1,  2,...,  Q. 


As  a corollary,  one  would  expect  that  subpopulations  with  greater 
susceptibility  to  the  disorder  should  be  tested  more  frequently.  This 
turns  out  to  be  the  case  provided  testing  costs  are  the  same. 

Corollary  3.1:  If  limD(t)  »°°  and  X^/c^  > X /c . , then  r*  > r*  . 

The  relation  provided  by  this  corollary  between  subpopulation  specific 
defect  incidence  rates  and  cost  per  test  on  one  hand,  and  the  relative  test 
ing  frequencies  for  the  various  subpopulations  on  the  other  hand,  does  not 
require  any  knowledge  of  the  shape  of  the  disutility  function  D(* ) , except 
that  it  is  increasing  to  infinity.  The  next  two  results  show  how  more 
specific  knowledge  of  D(» ) can  be  used  to  schedule  a mass  screening  pro- 
gram optimally.  The  first  provides  an  explicit  solution  to  the  program 
(4.5)  whetxD(t)  * atm  for  m > 0.  The  second  result  shows  how  the  solution 
for  the  case  D(t)  - at  serves  as  a bound  for  the  optimal  solution  when 
D(* ) is  any  convex  increasing  function. 


The  following  result  does  not  require  knowledge  of  the  exact  form  of 
the  function  D(* )--provided  that  it  is  differentiable  and  convex  increasing 

Theorem  4:  If  D' (• ) exists  and  is  increasing,  then  r*/r*  > / ^2C1^1C2 
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As  an  application  of  this  theorem  and  Corollary  3.1,  suppose  subpopu- 
lation #1  has  twice  the  incidence  rate  of  the  defect,  per  unit  population, 
as  does  subpopulation  #2.  Then,  if  testing  costs  are  the  same,  r*  < r*  < 
y/~T  r*  , provided  the  disutility  function  satisfies  the  hypotheses  of  the 
theorem. 

An  analogous  proof  to  that  of  Theorem  4 would  show  that  if  D' (• ) 
exists  and  is  decreasing  (hence  D(« ) is  concave  increasing),  then 

"ie  "fc 

r2^ri  < ^2C1^1C2  ^°r  ^]/ci  > ^2^C2’  ®ne  tec^nica^  difficulty  does 
arise,  however.  If  D'(»)  is  decreasing,  it  is  not  necessarily  true  that 
lim  D(t)  ■ and  hence  Theorem  3 cannot  be  utilized.  It  is  possible 

£ -•  CO 

to  get  around  this  problem  by  adopting  the  hypotheses  of  Theorem  3' . 

This  is  summarized  by  writing 

Theorem  4' : If  D' (•  ) exists  and  is  decreasing  and  if  the  optimal  solu- 

tion  (r^,...,  r ) to  (4.5)  exists  and  has  positive  components,  then 

r*/r*  < /T^T2  for  ^1^C1  > ^2/<c2  ' 

Kirch  and  Klein  obtain  essentially  this  result  for  c^  =■  c2  and  F a 

"ic 

uniform  distribution  over  an  interval  [0,  g]  . If  0 < r^  < g for  j * 

1,...,  Q,  then  their  F(*)  may  be  regarded  as  differentiable  and  increasing 
for  purposes  of  our  proofs. 

5 . Empirical  Estimation  of  D(» ) 

Although  data  may  not  be  available  currently  to  support  the  utiliza- 
tion of  particular  models,  that  should  not  deter  the  development  of  such 
models.  It  is  reasonable  to  expect  that  in  many  application  areas  in  the 
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future  much  additional  data  will  become  available;  for  example,  more  data 
will  become  available  concerning  the  stochastic  pattern  of  a particular 
disease's  development.  Furthermore,  a good  model  will  serve  as  a guide 
to  the  kinds  of  data  which  should  be  gathered. 

It  is  also  reasonable  to  expect  the  future  development  of  improved 
testing  technologies  capable  of  detecting  a disorder  in  a much  earlier 
stage  of  development  than  is  now  the  case.  Such  innovations  will  make 
screening  programs  more  attractive  and,  consequently,  make  more  impor- 
tant the  analytical  tools  to  design  such  programs  intelligently. 

However,  it  is  possible  to  provide  a simple  procedure  to  estimate 
empirically  the  disutility  function  D(* ) under  the  assumption  of  perfect 
detection.  Such  a procedure  is  necessary,  since  upon  the  detection  of  a 
disorder  there  may  be  no  way  of  directly  determining  how  long  the  disorder 
has  been  present;  although  that  length  of  time  could  not  exceed  x = ^ . 

It  is  assumed,  however,  that  at  each  detection  of  a defect,  the  degree  of 
disutility  incurred  due  to  that  defect  can  be  observed.  For  example,  in 
the  case  of  medical  screening,  at  the  detection  of  a tumor  its  degree  of 
development  can  be  noted  even  if  it  is  not  possible  to  determine  exactly 
when,  since  the  last  test,  the  tumor  originated. 

In  order  to  allocate  optimally  a screening  budget  among  differing 
subpopulations  (susceptibility  classes),  as  was  discussed  in  section  3, 
it  is  necessary  to  know  the  shape  of  the  disutility  function.  Under  the 
reasonable  assumption  that  D(* ) is  an  increasing  continuous  function,  the 
function  may  be  derived  in  the  following  manner.  For  a particular  popu- 
lation subject  to  a Poisson  defect  arrival  [S^}  with  rate  NX  , arbitrarily 
select  a value  for  x and  set  up  a prototype  mass  screening  program  for 


-25- 


the  population  in  which  tests  are  made  (using  a perfect  test)  at  times 
0,  x,  2x,...  The  number  of  defects  detected  at  the  time  jx,  which  are 
characterized  by  a disutility  less  than  or  equal  to  y,  record  and  desig- 
nate by  Q^(y).  Then, 

Qj (y)  is  observable  for  all  values  of  disutility  y > 0.  Further,  Q^(y) 
and  Qj(y)  are  independent  and  identically  distributed  random  variables 
for  i ^ j and  y > 0.  Therefore,  by  the  law  of  large  numbers, 

1 n 

lim  — £ Q.  (y)  = E[  Q (y ) ] for  all  y>  0. 

n-«  n J-l  J 

Hence,  the  function  y-*E[Q^(y)]  may  be  estimated  by  collecting  a 
sufficiently  large  sample  of  the  functions  , C^,... 

Let  H(y)  = EfQ^y)].  Once  H(»  ) is  estimated  in  this  manner,  the 
following  theorem  shows  how  the  desired  function  D(* ) may  be  obtained 
from  H(*  ) . 

Theorem  5:  If  D(0)  ■ 0 and  D(* ) is  continuous  and  increasing  on  [0,  x] , 

then 

D(t)  »■  H ^(NXt)  for  OS  t s x . 


6 . Summary 

In  section  3,  an  explicit  expression  for  the  expected  long-run  dis- 
utility per  unit  time  given  a screening  schedule  of  uniformly  spaced  test 
administrations  was  given.  A general  multi-population  screening  program 


was  also  formulated. 
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In  the  subsequent  sections  under  the  assumption  of  perfect  test  reli- 
ability, several  features  of  optimal  screening  schedules  for  a heteroge- 
neous population  have  been  developed.  Within  the  class  of  periodic  sched- 
ules, uniform  schedules  were  shown  optimal.  The  Kuhn-Tucker  conditions 
for  assigning  particular  optimal  uniform  schedules  (i.e.,  particular  test- 
ing frequencies)  to  the  respective  subpopulations  were  computed.  When 

lim  D(t)  » <=,  these  conditions  reduced  to  a simple  form  which  made  pos- 
t 00 


sible  the  following  more  specific  results:  (i)  conditions  (in  terms  of 

X and  c^)  under  which  one  subpopulation  should  be  tested  more  frequently 
than  another,  (ii)  explicit  testing  schedules  over  all  subpopulations  when 

^ A' 

D(t)  = at  , and  (iii)  bounds  on  r^/r^  when  D(* ) is  an  arbitrary  convex  in- 
creasing differentiable  function.  Finally,  a procedure  was  suggested  for 
the  determination  of  D(* ) under  perfect  test  reliability. 
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Appendix 

Lemma  A.  Let  {t  } be  a renewal  process  with  associa  d renewal  function 
G(* ).  And  let  (x!j)  , (Xk)  be  two  sequences  of  stochastic  processes  such 
that 

Ic  k v 

(a)  The  pair  of  processes  (X^,  X^)  is  stochastically  identical  to 
the  pair  (X* , X*)  for  k - 1,  2,... 

(b)  (X:j,  X^)  is  independent  of  Tk  for  any  j,  k. 

Then 


for  any  bounded  function  H. 

Proof: 

E £ H(Xk,  Xk  T k)  - £ EE[H(Xk  Xk,  Tk)|Tk] 
k-1  k-1 

- E EE[H(x),  X1  t k)  | t k]  . 
k-1  1 

The  last  equality  uses  hypotheses  (a)  and  (b).  Let 

f(s)  - E[H(xJ,  X*.  Tk)|Tk  - s] 

- E[H(x},  X*.  s)|rk  - s]  - E[H(xJ,  X*.  s)]  , 

where  the  last  equality  utilizes  (b).  Now,  since  H is  a bounded  function, 
so  Is  f.  Hence,  we  may  utilize  a result  in  finlar  (p.  286)  to  get 
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E S H(Xk  Xk  Tk) 
k-1 


2 Ef(TK) 
k-1 


CD  00 

f f (s)dG(s)  -J  E[H(xJ,  X*,  s)]dG(s) 


Proof  of  Theorem  1 : 


Q.E.D. 


m £ E[D<^,*-S%r(sk)l 

- 2^  e[  E[D(S^  A-Sk)|Sk]lj/r(Sk)] 


(A.l) 


Now,  since  ^ €{l/r,  2/r,...}, 

r >* 


E[D(Sk>A  - Sk)|Sk] 


- S D(—  - Sk)P(Sk  - n/r|Sk) 
n-1  r r,£ 


And  by  (3.3) 


P(Skx  - n/r|Sk) 


p(  ^ fYi(7'sk)  “ o}n(^<7-sk)  - i}) 


- p,<7-sk)  n [i-p,(±-sR)] 

It  lml  lr 


n-1 


ri  ok. 


k k 

by  (3.2)  and  the  independence  of  Y^(s)  and  Y^(t)  f°r  s f t. 
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Combining  Che  last  two  equations  gives 


E[D(S^£-Sk)|Sk] 


E D(^-sk)p^(S.sk)  n [i-p2(7-sk)l 

n«l  i-1 


This  equation  inserted  into  (A.l)  yields 


E[B  ] = Er  £ £ D(S-sk)p  (^-sk)  n [1-p/7-sk)]i1/r(sk)l 

r>Z,J  -k-1  n-1  r i-1  1 r j/r  J 


which  by  lemma  A 


J 2 D(S-»)Pj<a-s)  n <»)d» 


r ( j+l) / r * n-1 

NX  J S D(S- s)p  (*- s)  n [1  - p (i- s)]ds 

, . n-1  i-1 

j/r 


letting  u - — - s 


NX  E [ r D(u)p . (u)  n [l  - p.(u  - n ' 1)]du 

n-1  “ . i-1 

n-1  - 1 


letting  m - n - i 


® L n- 1 

E J r D(u)p^(u)  II  [ 1 - p^(u  - ^)]du 


n- 1-1 


which,  since  p(u)  - 0 for  u<  0 


■NX  E J r D(u)p^(u)  n [l  - p^(u  - ^)]du 
n"J+1  m“1 


i 
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The  last  equality  follows  since  for  m 2 v and  for  (v  - l)/r  £ u < v/r, 
m/r  a v/r  > u.  But  then  1 - p^(u  - rn/r)  - 1 for  m 2 V and  u within  the 
region  of  integration. 

Q.E.D. 


Proof  of  Lemma  1:  Since  D(* ) is  increasing 

ra  i*  Ya  rYa 

vj  D(s)ds  » j D(t/y )dt  > J D(t)dt  for  0 < y < 1 . 

0 0 0 

Now,  for  0 < a < 1 and  a * 1 - a , 

x y otx  j3  y 

a J D(t)dt  + g | D(t)dt  > J D(t)dt  + J D(t)dt  . 

0 0 0 0 

Q.E.D. 

Proof  of  Lemma  2:  For  an  arbitrarily  large  number  M,  it  is  necessary  to 

show  that  there  exists  a number  6 > 0 such  that  f(r)  i M for  all  r £ 6 . 

Since  lim  D(s)  - ®,  there  exists  a number  1/24  such  that  D(s)  i 2M 

s -*  ® 


for  s 2 1/2  5 . Then 
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1 l 

f (6  ) -6  J5  D(s)ds  i 6 J6  D(s)ds  i 6 (yg)2M  = M . 

o _L 

26 

Since  f(«)  is  a decreasing  function,  f(r)  a f(6)  «e  M for  r £ 5 . 

Q.E.D. 

Proof  of  Theorem  3:  Upon  taking  the  closure  of  the  feasible  set  in  (4.5) 

and  extending  the  objective  function  so  as  to  preserve  continuity,  the 
following  program  is  obtained: 

Minimize 
<rl»r2*-‘-*rQ) 

s . t. 

Since  the  objective  function  is  continuous  and  the  feasible  region  is 

“ic  "fc  ^ 

compact,  the  minimum  r * (r^,...,  r^)  exists.  By  the  strict  convexity 
of  the  objective  function  (Lemma  2),  r is  unique. 

The  linearity  of  the  constraints  and  the  differentiability  and  con- 
vexity of  the  objective  function  imply  that  the  Kuhn-Tucker  conditions 

•ff 

are  both  necessary  and  sufficient  conditions  for  r . The  Kuhn-Tucker 
conditions  are 

NjXJf'  (rJ>  +UlNjCj  'Uj+1  “ ° for  j - 1 Q 

and  (A. 2) 

uJ+1r*  - 0 for  j - 1, ....  Q . 

Mu  i 


£ N.\.f (r . ) 


J-l 


J j " J 


£ N c r s b 
j-l  J J J 

Tj  2 0 for  j - 1, ....  Q 
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Suppose 

r*  > 0 for  j - 1 Q . 


(A. 3) 


Then  uj+^  * 0 for  j ■ 1,...,  Q and  the  theorem  will  follow  from  the  first 
equation  of  (A. 2).  Of  course,  (A. 3)  is  also  necessary  for  feasibility 
with  respect  to  (4.5).  Hence,  it  only  remains  to  verify  (A. 3). 

If  for  some  i,  r*  • 0,  then 

Q * 

£ NX  f (r*)  a N X f(r  - 0)  - » 

J J J ill 

where  the  equality  follows  from  Lemma  2.  But,  it  is  easy  to  exhibit  a 
feasible  point  which  maps  to  a finite  value  of  the  objective  function;  for 


example,  (r. 


rQ 


) where  r^  =»  b 


b/?.  N‘c‘ 


for  j - 1,. 


This  con- 


tradicts the  optimality  of  r*,  hence  establishing  (A. 3)  and  the  theorem. 


Q.E.D. 


Proof  of  Corollary  3.1:  From  Theorem  3 

f ' (r*)/f ' (r*)  - (XjCj^)/ (Xj /Cj  ) - (Xj/cj)/(Xi/ci)  < 1 . 

Multiplying  each  side  by  fCr^),  which  is  a negative  number,  f'(r^)>  f'Cr^). 
The  corollary  now  follows  from  (4.7)  which  shows  f ' (• ) to  be  an  increasing 
function. 

Q.E.D. 


Proof  of  Corollary  3.2:  Since  lim  D(t) 

t -♦  Oft 

all  J - l Q 


XJf,(r?/Cj  " 


Theorem  3 provides  that  for 


cons  tant . 


(A. 4) 
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Now  compute 
f ' (r) 


L 


(*  r m 1 d T 

a 

-ma 

J atdtJ-di|_ 

n 

• (m+  l)rm  - 

. , . nrt-1 

(m+  l)r 

Therefore,  using  (A. 4) 

1 - [\Jf,(r*)/cJ]/[\1f’(r*)/c1] 


, *nrt-l 

. SVi 

* , *nri-l 

iVj 


j - 1 Q 


Hence, 


art-1 


’j 


rr  v J - 1 Q 

i j 


As  mentioned  above,  the  budget  constraint  in  (4.5)  is  binding;  therefore 

(nrt-1 ) /t — — - „ (nrt-1  V^-JT 

J-!  n J-l  7 X’ 


R * R 

b » £ N.c.r,  ■ r_,  S N,c, 


llcj  1 j-1  J V X1 


which  yields 


/q  ^yrsr 


Q.E.D. 


Proof  of  Theorem  4:  Since  D(* ) is  convex  increasing,  lim  D(t) 

0D 


Therefore,  by  Theorem  3 


X1f’(r*)/c1  - \2f'(r*)/c2 


f'(r*)/f(r*)  - X2C1/X1C2 


and  so 
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Hence,  to  prove  the  theorem,  it  suffices  to  show 

(r*/r*)2>  f ' (r*)/f ' (r*) 


that  is  (recall  f'  < 0) , 


*2  , . *. 
r-  f' (r.) 


r*2f,(r*) 


(A. 5) 


Now,  X^/^  > X2/c2  implies  by  Corollary  3.1  that  r*  > r*.  So  (A. 5),  and 

? 

thereby  the  theorem,  will  be  established  if  r f ' (r)  is  shown  to  be  an  in- 

d 2 

creasing  function;  i.e.,  — r f'(r)>  0. 


£r2f(r) 


^[r2(  JrD(s)ds* 

0 


l 

J r D(s)ds  - rD(~)^J 
0 


l 

2r  f rD(s)ds -D(-)  - D(-)+-D' (i) 
J r r r r 

0 


letting  x 


1 

r 


- - f D(s)ds  - 2D(x)  +xD'(x) 

X J 

o 


* — f xD(x)  - r sD'(s)ds"j  - 2D(x)+xD'(x) 

X |_  j J 

0 

2 r x 

* * x J sD' (s )ds + xD' (x) 

0 

> -£J  sD'(x)ds  + xD'(x) 

0 


since  x i s 

- -(2D'(x)/x)(x2/2)  + xD'(x)  - 0. 


Q.E.D. 
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Proof  of  Theorem  5:  Let  D(x)  ■ d.  Since  D(0)  - 0 and  D(* ) is  continuous 

and  increasing  on  [0,  x] , then  D ^ (•  ) exists  and  is  continuous  and  increas- 
ing on  [0,  d]  with  D *(0)  ■ 0 and  D *(d)  - x. 


by  Lemma  A 


Therefore, 


Et«iWl  - E[  £ l[o,y]CD<--sk>)lto,x] 

(Sk)] 

(*  X / \ 

* x / 

■"'J  1[0,y](D(x-t))dt  ' 

C 

o 

s. 

0 

0 

*NXJ  1[0,D-l(y)](s)ds  " Nx( 

x,  D_1(y)) 

0 

- NX  D_1(y)  for  OS  y £ d . 

D X(y)'  - H(y)/NX  for  y € [0,  d]  . 

Applying  D(« ) to  each  side, 

D^H(y)/Nx)  - y for  y€[0,  d] 


Let 


t ■ H(y)/N\  for  y€[0,  d] 


Then  t € [ 0,  x]  and 


y ■ H ^(NXt)  for  t€[0,  x] 


(A. 6) 


(A. 7) 


(A. 8) 


Inserting  (A. 7)  and  (A. 8)  into  (A. 6), 


D(t)  ■ y - H X(NXt  ) for  t € [0,  x] 


Q.E.D. 
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